Impact of grain evolution on the chemical structure of 

protoplanetary disks 

A.I. Vasyunin 

Max Planck Institute for Astronomy, Konigstuhl 17, D-69117 Heidelberg, Germany 

Department of Physics, The Ohio State University, 191 W. Woodruff Ave., Columbus OH, 

43210, USA 1 

OV vasyunin@mps.ohio-state.edu 

, — D.S. Wiebe 
<■ 

Institute of Astronomy of the Russian Academy of Sciences, Pyatnitskaya str. 48, 119017 
rn ' Moscow, Russia 

q ! dwiebe@inasan.ru 

22 ■ T. Birnstiel 

Max Planck Institute for Astronomy, Konigstuhl 17, D-69117 Heidelberg, Germany 
O ' birnstiel@mpia.de 

S. Zhukovska 

Max Planck Institute for Astronomy, Konigstuhl 17, D-69117 Heidelberg, Germany 

o: 

zhukovska@mpia . de 

> : 

^ , T. Henning 

Max Planck Institute for Astronomy, Konigstuhl 17, D-69117 Heidelberg, Germany 

henning@mpia . de 
CP. Dullemond 

Max Planck Institute for Astronomy, Konigstuhl 17, D-69117 Heidelberg, Germany 

dullemonOmpia . de 

ABSTRACT 



-2 - 



We study the impact of dust evolution in a protoplanetary disk around a 
T Tauri star on the disk chemical composition. For the first time we utilize 
a comprehensive model of dust evolution which includes growth, fragmentation 
and sedimentation. Specific attention is paid to the influence of grain evolution 
on the penetration of the UV field in the disk. A chemical model that includes 
a comprehensive set of gas phase and grain surface chemical reactions is used to 
simulate the chemical structure of the disk. 

The main effect of the grain evolution on the disk chemical composition comes 
from sedimentation, and, to a lesser degree, from the reduction of the total grain 
surface area. The net effect of grain growth is suppressed by the fragmentation 
process which maintains a population of small grains, dominating the total grain 
surface area. We consider three models of dust properties. In model GS both 
growth and sedimentation are taken into account. In models A5 and A4 all 
grains are assumed to have the same size (1CT 5 cm and 10~ 4 cm, respectively) 
with constant gas-to-dust mass ratio of 100. Like in previous studies, the "three- 
layer" pattern (cold midplane, warm molecular layer, and hot atmosphere) in 
the disk chemical structure is preserved in all models, but shifted closer to the 
midplane in models with increased grain size (GS and A4). Unlike other similar 
studies, we find that in models GS and A4 column densities of most gas-phase 
species are enhanced by 1-3 orders of magnitude relative to those in a model 
with pristine dust (A5), while column densities of their surface counterparts are 
decreased. We show that column densities of certain species, like C 2 H, HC 2 „+iN 
(n = — 3), H 2 and some other molecules, as well as the C2H2/HCN abundance 
ratio, which are accessible with Herschel and ALMA, can be used as observational 
tracers of early stages of the grain evolution process in protoplanetary disks. 

Subject headings: accretion disks, astrochemistry, opacity, ultraviolet: planetary 
systems 



Introduction 



The presence of numerous exoplanets and protoplanetary disks (PPD), as well as our 
own existence, strongly suggest that planet formation is ubiquit ous in the Milky Way (e.g., 
Mavor fc Quelozlll995l : iMavor fc Freill2003l : ludrv fc Santoll2007h . Still, it is one of the most 
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challenging problems for the modern astronomy to understand how planets form. This 
topic covers a tremendous range of micro- and macrophysics, and involves a wide variety of 
physical, dynamical, and chemical processes. 

According to the current paradigm, planets are assembled in protoplanetary disks sur- 
rounding yo ung stars, starting from the coag ulation of sub-micron dust grains into larger 



bodies (e.g.. iNatta et al.l 120071 ; iHenningl 120081 ). Circumstantial evidence indicates that this 



process is quite effective and rapid. For example, large-scale IR spectroscopic surveys of 
young stars in stellar clusters of various ages tightly constrain typical dispersal timescales of 

m ~ 1Q- 8 Mr, vr- 1 



(Haisch et al. 


2001 




Sicilia-Aguilar et a 


2006 


Sicilia-Aguilar et al. 


2009: 


Fedele et al. 


2010) 



Hernandez et al. 2008; Oliveira et al 



2009 



Inner, planet-forming zones of some older 
disks are cleared of sub-micron-sized dust grains within just ~ 1 Myr (e.g., I Graham et al. 
2007l : lDutrev et al-lbooa ISalvk et aDl2009i Ehalmann et al.lboioh . Furthermore, the analysis 
of the mineralogical and chemical composition of unaltered chondritic meteorites shows that 
in the Solar Nebula different types of meteorit es have condensed and formed via grain agglom 
eration within a period o f 
19941 ; iThrane et alJliooS ). 



2 — 5 Myr (e.g.. IWasson fc Kallemeynlll988t iPodosek fc Cassen 



However, observational data on grain growth are still scarce and limited to the over- 
all spectral energy distribution (SED) parameters. The prese nce of micron-sized gr a ins in 



disks has been inferred from mid- infrar e d spectroscopy (e.g., iBouwman et al.l l200ll 12008 



Apai et al.l l2004t Ivan Boekel et al.l l2004t iKessler-Silacci et al.l 120061 ). A number of nearby 
protoplanetary disks has been imaged with the VLA and Australia Telescope Compact Ar- 
ray (ATCA) at millimeter and centimet er wavelengths, showing evidence of significant grain 



grow t h up to at least pebble (cm) sizes dTesti et al. 



2006; Natta et al. 2007; Cortes et al. 2009 



2003 



Lommen et al. 



Wilner et al. 2005; Rodmann et al. 



20091 ). On the other hand, there 



can be other observables indicative of the grain ensemble evolution. 

In essence, grain evolution is controlled by collisions, leading to coagulation and frag- 
mentation of dust particles. For mm-sized particles the bouncing effec t comes into play , 
making the process more complicated and further growth less efficient (jZsom et al. 2010 ). 



Radial drift and sedimentation to the disk midplane become important if the growth of grain 
particles proceeds further than about 1 cm. Grains evolve faster in the dense disk interior 



unchanged (ISM-like) for the entire life of the system ( 


Schmitt et al. 


199' 


I; Weidenschilline 


1997; 


Ciesla 


2007; 


Ormel et al. 


2007; 


Brauer et al. 


2008a; 


Birnstiel et al. 


2009. 


2010 


). It is 



therefore reasonable to expect that the evolution of grain sizes and spatial distribution is 
somehow reflected in the physical and chemical structure of the disk. 



-4- 



In current models of disk structure it is assumed that the temperature distribution 
within the disk, and thus its structure, is controlled by the radiative transfer of energy 
generated due to viscous dissipation and absorption of the stellar radiation. In the disk 
atmo sphere, photoelectric heat i ng from dust grain s becomes important for energy balance 



[e.g., IGorti fe riollenbachl 120091 ; IWoitke et ail 120091 ). Dust opacities, which are changed by 
grain growth and sedimentation, are the crucial element in this picture. Also, because of dust 
evolution, stellar UV photons penetrate more easily into the disk interior, heating it, and 
thus preventing rapid freeze-out of molecules and enhancing photodestru ction of gas-phase 
chemical species (e.g.. lAikawa &: Nomurall2006t iJonkheid et al.ll2006l 120071 ). Local variations 
in the dust-to-gas ratio, dust density and average grain size may all affect the efficiency of 
gas-grain interactions, and thus the efficiency of freeze-out and surface chemistry. 

The problem of possible interrelations b etween dust evo l ution and chemistry has been 
addressed in the literature a number of times. IJonkheid et al.l (120041 ) studied the importance 
of dust settling on the disk chemical evolution, assuming that the dust size distribution 
remains the same along the disk radius. It was found that dust sedimentation has a great 
impact on the gas temperature in the disk and that dust and gas temperatures are signif- 
icantly different in the optically-thin disk region. On the other hand, sedimentation was 
found to not strongly influence the gas-phase chemistry in the model considered due to the 
assumption that PAHs are still well-mixed with the gas and therefore efficiently absorb stel- 
lar and interstellar UV photons. At the same time, elevated temperatures of the upper disk 
layer in the model with grain sedimentation lead to excitation of higher-lying transitions. 



Aikawa fc Nomura (120061 ) considered the effect of dust growth on chemical abundances 



in a well-mixed disk around a low-mass star. They modeled dust and gas temperatures in the 
disk self-consistently and adopted a power-law grain size distribution. The grain growth was 
simulated by raising the maximum grain size up to 10/xm. An extended chemical network 
with gas-grain interactions and surface processes was applied. They found that the main 
effect of dust growth is a shift of the molecular layer towards the disk midplane. At the 
same time, column densities for many species were not affected by grain growth. Finally, 
they identified a few sensitive tracers of grain growth in disks, such as HCO + , ILt, H 2 D + . 



Later, IJonkheid et al.l (120071 ) modeled chemistry and gas temperature of evolving Herbig 
Ae disks with decreasing mass or dust settling. They utilized a 2D radiative transfer code 
to model gas/dust temperature and photodissociation rates. It was found that the disk 
chemical structure strongly correlates with the disk mass: while some photofragile molecules 
like HCN are abundant only in massive disks, their fragments, like CN and CCH, become 
abundant in low-mass disks. Furthermore, dust settling affects the gas temperature in the 
disk atmosphere, and consequently abundances of disk atmosphere tracers like O and C + . 
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They identified that line intensity ratios of CO/ 13 CO, CO/HCO+ and [01] 63/im/146/im can 
be used to distinguish between disks undergoing grain settling and photoevaporated disks. 

In the present work, we move one step further and omit the simplifying assumptions 
that grain evolution proceeds uniformly over the entire disk, or that it can be described by 
artificially varying the grain number density and upper grain size limit. The goal is to isolate 
possible molecular tracers of grain evolution in the situation when grain size distribution and 
dust-to-gas mass ratio vary smoothly in the disk, as predicted by a dust evolution model. The 
model used includes coagulation, fragmentation, and sedimentation , assuming equilibrium 



betwe en turbulent stirring and gravitational settling, as proposed in iDullemond fc Dominik 



(120041 ) . No bouncing effects are included. We focus on the combined effect of dust growth 
and settling on the chemistry in a protoplanetary disk around a low mass T Tauri star. A 
large gas-grain chemical network with surface reactions is utilized. This is the first paper in 
a series, where we intend to study various aspects of interrelations between grain evolution 
and chemistry. 

The outline of the paper is the following. In section [2J the setup we constructed to 
model the physical structure of disk with evolving dust is described. Section [3] contains 
a description of the detailed chemical model used to simulate the chemical structure of 
the protoplanetary disk. In section H] we present modeling results. In paragraph 14.11 dust 
evolution due to growth, fragmentation and sedimentation is discussed briefly. In paragraph 
14.21 change to the UV field in the disk caused by evolution of dust properties are discussed. 
Paragraph 14.31 is devoted to the analysis of the impact of grain evolution on the chemistry 
of a protoplanetary disk. In section HJ we discuss whether or not the chemical composition 
of a protoplanetary disk can be used as a tracer of grain evolution, and in which ways our 
model can be improved. Section [6] contains the conclusions of this study. 



2. Description of the model 

A protoplanetary disk, even in a simplified description is a complicated evolving system 
controlled by a large number of interrelated physical and chemical processes. However, given 
the uncertainties and deficiencies in the physical description of such a system, it would be 
impossible and impractical to try to build a self-consistent dynamical model coupled to grain 
evolution and disk chemistry. For this reason, we have selected an approach which is by no 
means self-consistent, but is relatively straightforward to implement and allows us to isolate 
the purely chemical effect of the dust growth. 

The disk chemical modeling in this study comprises of four independent steps. First, we 
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need a template of a disk physical structure, i.e. gas density and temperature distributions. 
These distributions are obtained from a 1+1D steady-state accretion disk model. In the 
second step, the resulting thermal and density structure is used as an input for the grain 
growth model. The third step is to simulate the evolution of the dust grain ensemble in 
the disk midpl a ne du e to gr ain coagulation, fragmentation and sedimentation according to 



Birnstiel et al.l (120091 120101 ). Third, the vertica l distribution of dust i s recon structed in 



stirring-mixing equilibrium based on the work by iDullemond fc Dominikl (120041 ) . 

Finally, the grain parameters (size distribution at each disk location, as well as dust-to- 
gas ratios) after 2 Myr of evolution are used in the simulation of the disk chemical structure. 
These steps are described in more details in the following subsections. 



2.1. Physical structure of the disk 



A popular formalism, which is used in many studie s of di sks surrounding young stars, 
is based on the so-called a-model of IShakura fc Sunyaevi (119731 ). In this model it is assumed 
that the transport of angular momentum in the disk is caused by turbulent viscosity, with 
the viscosity coefficient v expressed as aHc s , where H is the disk semi-thickness at a given 
radius, c s is the isothermal sound speed, and a is the viscosity parameter, which is usually 
assumed to be of the order of 10 -2 in protoplanetary disks. 

As a first approximation for calculating the disk thermal structure of the disk, we assume 
thermal coupling between dust and gas. This assumption is well justifie d for optical depths 
r > 1, as was shown by t he disk models with detail thermal balance (lAikawa fc Nomura 
20061 ; IJonkheid et al.ll2006l ). Sedimentati on and coagulation o f grains influence the disk tem- 
perature via a decrease in dust opacities (IBirnstiel et al.ll2009l ). which makes the atmosphere 
more transparent and hotter, and shifts the r = 1 surface deeper into the disk. The investi- 
gation of the response of the disk thermal structure on dust evolution is a separate problem, 
which will be addressed by future papers in this series. 

The disk structur e is c alculated with a simpl ified version of the disk model described by 
D'Alessio et al.l (119981 ) and lD'Alessio et al.l (11999I ). Unlike D'Alessio et al., in the calculation 
of gas temperature we take into account only two heating sources, namely, viscous dissipation 
and disk irradiation by the central star. Comparison of our disk models with the D'Alessio et 
al. results shows that this simplification does not significantly alter the derived disk physical 
structure, and thus modeling of the disk chemistry is not affected. 
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The following equations are used to describe the disk vertical structure: 

dP 

dT _3npF_ 

dz 4acT 3 ' 1 ' 

dF 9 

— = -aPn(R) + T ilI . (3) 

Here P is the gas pressure, g is the acceleration due to gravity, p is the gas density, T is the 
gas temperature, F is the flux of the thermal disk emission, Q(R) is the disk (Keplerian) 
angular velocity at a distance R from the central object. The first term on the right hand 
side of Eq. ([3]) gives the viscous heating of the disk, while r irr is the heating rate due to disk 
irradiation by the central object. 

The two-dimensional structure of the disk is modeled in the so-called 1+1D approxima- 
tion. It is assumed that at each radius the disk is in hydrostatic equilibrium (Eq. [T]). Thus, 
the vertical structure of the disk can be found from the solution of equations ([IH3]) indepen- 
dently at each R. To calculate the disk temperature at each vertical slice, we integrate the 
transfer equation ©. The disk itself is described by the following parameters: 



• accretion rate M; 

• viscous parameter a; 

• temperature, radius, and mass of the central object i?*, M* (radiation is assumed 
to be that of a black body); 

• inner and outer disk radii i?; n and -R ut- 



The outputs of the model are gas surface density distribution S ffas (i2), and the 2D 
density and temperature structure of the gas. The strengths of the UV and X-ray fields over 
the disk can then be calculated using the information about the gas and dust properties in 
the disk. 

Using this model, we simulated a DM Tau-like system. The system has a disk with an 
outer radius of 550 AU, an inner radius of 0.03 AU, an accretion rate M — 2 ■ 10~ 9 M Q yr _1 , 
and a viscosity parameter a = 0.01. The disk is illuminated by UV radiation from the 
central star w ith an intensity x = 500 at R = 100 AU in units of the mean UV field of 



Draind (119781 ) and by the interstellar UV radiation. Two representations for the stellar UV 



field are used. One is the scaled Draine field, and another one is the representative observed 
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spectrum as discussed by iBergin et al.l (120031 ) . The effect of these different radiation fields 
will be discussed in section 3.1. 



The star is assumed to have an X-ray luminosity of 10 30 erg s _1 (IGlassgold et al.lll997l ). 
X-rays induce the same reactions as cosmic ray particles, so that X-ray ionization rate is in 
fact simply added to cosmic ray ionization rate. The disk mass is M = 0.055 M . The disk 
surface density profile and its power-law approximation are depicted in Fig. [U The best-fit 
line has a sl ope of 0.84, which is somewhere in b etween of the Minimum Mass Solar Nebula 
index (1.5) (1 Weidenschillingl 1 19 77t lHavashilll981l) and the value derived for PPDs from SED 
at millimeter wavelengths, (> 0.5) ([Andrews fc Williams! 120071 ). The density and thermal 
structure of the disk are shown in Fig. [2j The main model parameters are listed in Table [TJ 
We do not consider a disk region beyond i? ut- 



2.2. Dust evolution in the disk 



Sever al authors have constructed theoretical models for grain growth in protoplanetary 



disks (e.g. iNakagawa et al.lll98ll; ISchmitt et a 



2005: Qrmel et al. 2007: IZsom k Dullemond 



. 1997 



2008 



Weidenschil 



Brauer et al. 



mg 



1997; Dullemond fc Dominik 



2008bl ). The main chal- 
lenges of grain growth modeling are the following: first, growth from ~ 0.01 fim grains to 
centimeter-sized grains encompasses many orders of magnitude. A growth simulation must 
be able to conserve mass, since large particles can grow by accreting large numbers of very 
small grains. Second, particle size is not the only parameter determining the properties of 
grains: porosity, composi tion and the various possible outcomes of collisions further extend 



the parameter space (see ISchmitt et al.lll997l : iGiittler et al.l 120 101 ) . 



Her e, we use a sli g htly m odified version of the code presented in iBirnstiel et al.l (120101 ) 
(see also iBrauer et al.l ( l2008bh ). This is a statistical, mass-conserving code which implicitly 
solves the Smoluchowski equation, taking coagulation, fragmentation and cratering into ac- 
count. It is important to note that we ignore radial drift of dust particles in the present 
study. 

The mathematical formalism of the model is the following. The number density distri- 
bution n(m,r,z) is a function of mass m, radius r, height z above the midplane, and time. 
We do not consider grain porosity as an additional parameter. We define the vertically 
integrated number density per mass m as 



N(m, r) 



n(m, r, z) dz. 



(4) 



We assume that coagulation and fragmentation Kernels K and L (defined below) are 
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independent of z. We can now describe the time-evolution of this distribution by a vertically 
integrated version of the Smoluchowski equation, 



dN{m) 

at 



+ Jo" 1 ^ 2 N(mf) N{m — m') K(m', m — m') dm' 
— J °° N(mf) N{m) K(m, m!) dm 1 



(5) 



+ | J Jq 00 N{m') N{m") L(m', m") S(m, m', m") dm' dm" 

— J °° N{m') N(m) L{m, m') dm' , 

where the radial dependence has been omitted, since we treat each radius independently, 
neglecting radial movement of dust. The right-hand side terms of Eq. [5] (from top to bottom) 
correspond to gain and loss by coagulation and gain and loss through fragmentation. 



The coagulation kernel K(m, m') and the fragmentation kernel L(m, m') are then given 



by, 



K(m, m') 
L(m, m') 



^2Ti{h(m) 2 + h{m') 2 ) 
1 

^/2ir(h{m) 2 + h{m') 2 ) 



p c ■ a(m, m') ■ Au(m, m') 
Pf ■ a(m, m') ■ Au(m, m'), 



(6) 



where p c and pi are the coagulation and fragmentation probabilities, respectively, which have 
a sum of unity, Au(m,m') is the relative particle velocity, and a(m,m') is the sum of their 
geometrical cross se ctions. In this work, we co nsider Brownian motio n, vertical settling (see 
Brauer et al.ll2008bl ). and turbulent motion (see lOrmel fc Cuzzill2007l ) as physical effects that 
produce the relative particle velocities. 

Particles colliding with a relative velocity higher than the critical velocity itf are assumed 
to fragment into a pow er-law size distribution of fragments (i.e., S(m,m' ,m") oc m~ L83 , 



sec 



Brauer et al.l l2008bl ) if the particle masses differ by less than one order of magnitude. 



Otherwise, the smaller body is assumed to excavate mass from the larger one by cratering, 
where the amount of excavated m ass equals the mass of the smaller body. The fragmentation 
velocity ui is taken to be 10 m/s (IBirnstiel et al.ll2009l ). 



This model is related to a model of gas disk structure through two input parameters: 
radial midplane temperature distribution T mi( j(-R) and the radial gas surface density distribu- 
tion S gas (i?). A more comprehensive description of the physics of coagula tion /fragmentation 



and of its numerical implementation can be found in lBrauer et al.l (j2008bl ) and IBirnstiel et al. 

(boioh . 
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The numerical model described above provides us with surface densities of grains of 
different masses. To calculate the vertical distribution pi(R,z) of grains, we assume that 
it is controll ed by the equilibrium between gravitational settling and turbulent stirring, as 
proposed by iDullemond Dominikl (12004 ) . Vertical scale heights are computed separately 
for each size bin. 



2.3. Average grain size for the chemical model 

An output of the dust evolution model is the grain size distribution at each disk location. 
In principle, this distribution can be incorporated into the chemical model by considering 
surface reactions and gas-dust interactions for each size bin separately. However this would 
make the model too complicated, and we elect to use a single grain size a(r, z) in the chemical 
model, which is computed from the equilibrium local grain size distribution in each disk 
location. The averaging is performed in a way which preserves the total surface area and 
the total dust mass. For a grain size distribution f(a, r, z), where f(a, r, z)da is the number 
of grains with radii between a and a + da in the disk point with circumstellar radius r and 
height above midplane z, this means that 



, ma x 



4 . 4 i 

-%a 3 (r, z)n = -n j f(a,r,z)a 3 da, (7) 



47ra 2 (r, z)n = An / f(a,r,z)a 2 da, (8) 

where n is the total number density of "representative" equal-sized particles. By dividing 
Eq. (J7J) by Eq. (jSJ) one obtains 

f(a,r,z)ada 

a{r,z)= ;-J) 2 ■ (9) 

I j(a,r, z)a z da 

An ensemble of grains with the mean radius a(r, z) has the same mass and total surface 
area as the original ensemble. We use a(r, z) exclusively in our astrochemical model described 
below when dealing with gas-grain interactions and surface chemical processes. 

It is important to note that the way we calcula te the averag e grain s ize is different 



from t hose usually used in works on grain growth (e.g.. iBrauer et al.l (l2008bl ); iBirnstiel et al. 



(120101 ) etc.). In such works, the averaging is mass-weighted, i.e., emphasis is put on the size 
of a small fraction of the most massive grains. In this paper, on the other hand, averaging 
is surface-weighted, i.e. emphasis is put on a population of small grains which dominate the 
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grain surface area as well as the total number of grains. In other words, if simulations of grain 
growth would give us a distribution of grains in which a small number of big boulders co- 
exists with a big number of tiny grains, the average grain size of such a distribution calculated 
with expression (jH]) will be close to the size of small grains, contrary to the "usual" definition 
of average grain size. 



Chemical model of the disk 



In principle, simultaneous modeling of the disk chemical evolution, together with grain 
evolution, requires a non-local approach. G rains not only grow , but also move within the 



disk due to sedimentation and radial drift (IBrauer et al.l l2008bl . e.g.), bringing with them 



molecules in icy mantles. However, a coupled chemo-dynamical disk model with grain growth 
would be computationally expensive, and difficult to analyze from the chemical point of view, 
so in the present study we choose to neglect transport effects. Thus, we model a chemical 
structure in the disk, taking dust properties and physical conditions from the previously- 
described disk model, and assuming that these conditions do not change over the timespan 
of chemical evolution considered. This timespan is taken to be 2 Myr, but the particular 
choice does not significantly affect our conclusions, as the abundances of most species are 
equilibrated very fast. 



In this study, we utilize the chemical model d escribed in IVasyunin et al.l (120081 ) . Initial 
elemental abundances are taken from Table 1 in IVasyunin et al.l (120081 ). Compared with 
the original chemical model, we made several modificat ions. First, dissociative recom bina- 
tion of ions on grain particles is included according to lUmebayashi fc Nakanol (119801 ) and 
Semenov et al.l ( 120041 ) . Second, for gas-grain chemical interactions and grain surface reac- 
tions, we used variable dust-to-gas mass ratios and the average grain size calculated at each 
disk point usin g data from the grain growth model. For grain surface reactions, we utilized 
model H from IVasyunin et al.l (120091 ). In this model, only thermal hopping is a source of 



mobility of surface species, and a high diffusion-to-desorption energy ratio of 0.77 is adopted 
for all species. Under these conditions, stochastic effect s in grain surface chemistry are neg- 



ligibl e and classical rate equations may be safely used (IVasyunin et al.l l2009t iGarrod et al. 



20091 ). In our model molecular hydrogen is formed via the surface reaction H + H — > H2. 
The rate of this reactio n, as well as of other su rface reactions, is calculated according to the 
formalism described by lHasegawa et al.l ( 119921 ). The rates are inversely proportional to the 
surface area of an individual dust grain. 



The cosmic ray ionization rate has been adopted according to Eq. 19 from ISano et al. 



( 120001 ). The unattenuated cosmic ray ionization rate was taken to be 1.3 • 10 



- 17 s- 1 . 



The 
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final and most important change is the detailed treatment of photoprocesses which take into 
account the shape of the UV spectrum of the central star, and its attenuation in the disk 
due to absorption of the UV photons by dust grains (see Section [XT]) . Below we discuss in 
detail the treatment of these photoprocesses. 



3.1. Photoreactions 



Photoreactions represent the crucial element of chemical networks designed for the upper 
disk atmosphere, and it is this element that is most affected by the dust evolution. In the 
conventional approach, photoreaction rates are estimated as 



k ph = x^oexp(-7^ 



v 



(10) 



where k is the reaction rate for the unshielded interstellar UV radiation field, x is the 
UV intensity scaling factor, and 7 is a parameter used to account for the different dust 
extinction in the UV and visual wavelengths. Generally speaking, Ay estimated from the 
gas column density cannot be used in this case, because it is based on a particular choice 
of dust properties (opacities and dust-to-gas mass ratio), and on the assumption that these 
properties do not vary in the medium. In our study dust, parameters are position-dependent, 
so we have to be more careful. To compute rates of photodissociation and photoionization 
one uses the general expression: 



k ph (r, z) 



An 



Iy{r, z) 
hv 



o v dv, 



where the mean UV field intensity I v {r, z) is obtained as a solution of the radiation transfer 
problem. In this study we adopt a simplified approach, in which I u (r,z) is given by the 
equation 

I v {r,z) = i: t (r,z)ex P (-T:\r,z))+I»exp(-Tl s (r,z)), (12) 

where I^(r,z) is the unatten uated diluted stellar radiation field, I® is the so called Draine 
radiation field (jDraindll978l ). We investigated two representations for the stellar radiation 
field. In the first representation the observed spectrum for TW Hya is used, scaled in such 
a way that the integrated UV int ensity is about 500 in units of the Draine field at 100 AU 
from the sta r (jBergin et al.ll2003l ). A combination of FUSE (IHerczeg et al.ll2002l ) and HST 
(lAyresll2010l ) data is utilized, and the resul ting spectrum is sho wn in Fig. [SJ Essentially, this 
spectrum is very close to the one used in iBergin et al.l (120031 ) . In another representation, 
the stellar radiation is assumed to be given by the scaled Draine interstellar field, with the 
same normalization. Comparison of results for the two representations shows that both 
spectra lead to the same general conclusions about the chemical structure of the disk both 
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with e volved and unevolved dust. This is consistent with findings of Ivan Zadelhoff et al. 



( 120031 ). These authors compared molecular distributions and column densities for the cases 
when the UV spectrum is given either by the scaled ISFR spectrum or the smoothed TW 
Hya spectrum and found no significant differences (their models A and B). In our case, the 
relative similarity of results for the two representations is caused by effective attenuation of 
the stellar radiation in the disk, so that the radiation field is in both cases dominated by 
the interstellar field everywhere in the disk, except for the uppermost layers. The situation 
would have been different if scattering were present in our model. When the scaled ISRF 
is used to represent the stellar radiation, scattering is in some sense equivalent to just a 
somewhat different scaling factor for the vertically penetrating radiation. However, when a 
more realistic stellar spectrum is used, like the La-dominated TW Hya spectrum, scattered 
stellar radiation would n oticeable change p hotoreactio n rates for species w i th lar ge cross- 



van 



Zadelhoff et al. r|2003h did find 



sections in the La range (jBergin et al.ll2003l ). Note that 
a significant dependence on the spectrum of the radiation field in a more general sense, that 
is, for the case when no UV excess is taken into account at all (their model C). 



Also, Ivan Zadelhoff et al.l ( 120031 ) investigated how different approaches to the radiation 
transfer, namely, full 2D and simplified calculations, such as described below, impact the 
chemical composition of a protoplanetary disk. They found that column densities of most 
species are not sensitive to the chosen approach. Important exceptions are the molecules 
CN, C2H, CS and HCN, but only in the inner disk. Their column densities do significantly 
depend on the adopted UV radiation transfer model. However, in our model, the strength 
of the stellar UV field is m ore than an order of magnitude weaker than that considered by 
van Zadelhoff et al.l toO$ ). 



The optical depth r^(r, z) is computed along the ray between the star and the current 
location, so that 



_st 



n u (r, z)p(r, z)ds. 



(13) 



Dust sedimentation is taken into account in the dust density p(r,z). The absorption coeffi- 
cient K l ,(r, z) per unit dust mass is calculated at each location along the ray as 



K u (r, z) — C J daf(a,r, z)Q(a,b , )iia 2 , 
where the normalization coefficient C is defined by the condition 

1 = C J daf(a,r,z)^7ca 3 p d . 



(14) 



(15) 



Absor ption efficiency factors Q(a, u) are taken from lDraine &: Led (119841 ) and lWeingartner &: Draine 
( 1200 ll ). For the purpose of opacity computation, we assume that dust consists of silicate and 
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carbon (graph ite) components contributing , correspondingly, 70% and 30% to the total dust 
mass density (jWeingartner &: Draind 1200 ll ). In Eq. [15] pd is the density of dust material, 
taken to be 3 g cm -3 for silicate grains and 2.24 g cm -3 for graphite grains. In the expres- 
sions above, f(a,r,z) is the local evolved grain size distribution that comes from the dust 
growth and sedimentation model. Considering coagulation and fragmentation, we do not 
distinguish between the two materials as experimental data show a similar dust evolution 
behavior for both types of grains (J. Blum, private communication). 

The optical depth r* s (r, z) is computed in the same way, but in vertical direction from 
the disk surface to the height z. 



The reaction cross sections are taken from Ivan Dishoeck et al.1 (120061 ) and supplemented 
with the additional data for the Lya wavelength, with the exception of the photodissociation 
cross sections for NO, HCN, N0 2 , S0 2 , OCS, NH 3 , CH 4 , H 2 2 , and C 2 H 2 . For these 
molecules, monochromatic photodissociation cross sections from the AMOFtj database have 
been used in our study. 

To compute the rates of photoreactions for which wavelength-dependent cross-sections 
are not available, we use the standard expression fflOl) . however with the modified procedure 
for computing Ay. The conventional expression for optical extinction 



Ay = JVh/1.59 ■ 10 21 cm 



(16) 



is valid for dust which is well-mixed with gas and consists of single-sized grains (10~ 5 cm) 
and a dust-to-gas mass ratio 0.01, but it can easily be re-normalized for other values of these 
parameters. We integrate Ay along the ray between the star and the considered location (or 
in a vertical direction for interstellar radiation), taking into account the variable dust-to-gas 
ratio along the ray, and a variable average grain size a(r, z). 



Self-shielding for H 2 photodissociation is computed using the iDraine fc Bertoldil (119961 ) 
formalism, with the modified Ay value used to account for dust attenuat ion. The s e lf- an d 
mutual shielding for CO photodissociation is computed with data from iLee et al.1 ( 119961 ). 
Their dust attenuation is not included in the resultant shielding value, as it is replaced by our 
optical depths. To compute self-shielding factors, one would need to know the CO column 
density along the ray, which in the chemical study would entail iterations. We adopted 
a much simpler approach, assuming that CO self-shielding is proportional to the overall 
column density with an average abundance of 10 -6 . To check whether this simplification 
affects our conclusion, we run two models, bracketing possible self- shielding values. In one of 
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the models, maximum self-shielding is assumed which would correspond to all C atoms being 
in CO molecules. In the other model, CO self-shielding is assumed to be zero. Comparison 
of these two models showed that various treatments of self-shielding do not change our 
main conclusions while they do change some numbers, so we believe that the use of this 
simplification is acceptable. Anyway, species that turned out to be sensitive to the CO 
self- shielding treatment (CS and H 2 CO) are not included in Tables |2] and [3j 

In the model, several important simplifications are made. First, we assume that dust 
and gas temperatures are equal everywhere in the disk. Next, we neglect transport processes, 
such as radial and vertical mixing, and their influence on molecular abundances. Finally, 
when simulating grain-surface chemistry, we utilize the single average grain size instead of 
the real grain size distribution. Impact of these simplifying assumptions on the results of 
this study is discussed in Section 



4. Results 
4.1. Dust evolution in the disk 



Using the model of grain growth and fragmentation described in Section I2.2[ we calcu- 
lated the evolution of surface densities f or grains of differen t sizes, starting from the standard 



MRN grain size distribution at t — 0. ( IMathis et al.lll977l ). Snapshots of surface density as 



a function of grain size for four time moments are shown in Fig. [3J The corresponding verti- 
cally integrated average grain size calculated according to Eq. using vertically integrated 
number densities from Eq. @ is shown in Fig. HI 

The two peaks in Fig. H] appear because small grains typically grow together in a more 
or less uniform bump moving to larger sizes. So, the whole distribution shifts up in size 
until grains reach the fragmentation barrier. Once that happens, most of them fragment, 
and the small sizes get re-populated. The average grain size drops again, and an equilibrium 
between growth and fragmentation is established. 

Since the disk density is higher closer to the star, grain growth proceeds faster there (in 
our model we do not consider the radial drift of grains). One can see that at R < 30 AU 
a steady-state in the grain surface density distribution is reached within less than 10 4 yrs. 
However, in the outer disk (R > 100 AU) grain surface density evolves much slower. A 
steady-state distribution is barely reached there within 10 6 yrs. In order to ensure that 
grain parameters are equilibrated at all radii, we ran the grain evolution model for 2 Myr. 



It is important to point out that the highest surface density (Fig. [3]) corresponds to 
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grains which are bigger than the average particle size, as shown in Fig. HJ It is in these 
large grains that most of the dust mass is contained, while small grains dominate the dust 
surface area. Therefore, the presence of small grains in the disk mediates the net effect of 
grain growth from the point of view of gas-grain interactions, as they depend primarily on 
the available grain surface. 

In general, the grain growth is quite modest in the modeled disk. The largest increase in 
the surface averaged grain size is achieved in the inner disk midplane (10 < R < 100 AU) 
with a mean grain size of 10~ 4 cm (1 /iin). In the outer disk region (R > 100 AU) the 
average grain radius increases relative to the initial value by a factor of only 3-5. 

In contrast, differential dust settling seems to be much more significant. 2D distributions 
of the total dust-to-gas mass ratio and average grain size are shown in Fig. Dust settling 
reduces the dust-to-gas mass ratio in disk regions above the midplane. Everywhere above 
Z/R ~ 0.1 this ratio is significantly below the canonical value of 0.01. In the chemically rich 
intermediate layer it varies between 0.001 and 0.01, while in the disk atmosphere it drops 
down to 10 -4 — 10 -5 . At the same time, in the midplane the dust-to-gas mass ratio can reach 
a value of 0.02. Similarly, the average grain radius decreases above the midplane, though 
its change is not that dramatic. In the disk atmosphere, the average grain radius is only 
10 -6 cm (0.01 /im), because large grains efficiently settle down to the midplane. 

Summarizing both the vertically- averaged and 2D parameters of the evolved dust, we 
note that the overall effect of the grain growth and sedimentation in our disk model is the 
global decrease of the total surface area of dust grains (see Fig. [6]). This is less pronounced 
in the outer part of the disk (R > 100 AU), corresponding to a factor of only 2-3. In the 
inner disk region (R < 100 AU) the decrease of the dust surface area is more prominent, 
i.e., up to one order of magnitude. The ratio of evolved to original grain surface area has 
a steeper radial dependence than variations in the vertical direction. This is because the 
decrease in total dust surface area above the midplane, caused by the dust settling, is partly 
compensated by the presence of small particles dominating the surface area. 

4.2. UV field in the disk 

For the radiation field in the disk, we consider three basic models. In model A5, all 
grains are assumed to have a single size of 10~ 5 cm. This is the model of choice in most 
astrochemical computations. Model A4, in which all grains have the same radius of 10 -4 cm, 
can be considered a rough representation of grain growth (without sedimentation). Finally, 
in model GS both grain growth and sedimentation are taken into account fully. For brevity 
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we will call the GS and A4 models "evolved" 
"standard" model. 



while the A5 model is referred to as the 



In all models, the disk is illuminated by both stellar and interstellar UV radiation. For 
disk chemistry, the strength and shape of the UV field is of primary importance, as it controls 
the rates of photoionization a nd photodissociation, as well as the efficiency of the photodes 



orption of grain mantles (e.g., iBourdon et al.lll982t IWestley et al.lll995l : Ivan Zadelhoff et al. 



20031 ; lOberg et all 120071 ). Due to the presence of dust, the UV field is attenuated as it pen- 
etrates deep into the disk. As mentioned in the introduction, the attenuation of starlight 
in the disk atmosphere can also be caused by PAH particles. IR observations of protoplan- 
etary disks show that PAHs are an abundant component of some disks around Herbig Ae 
stars 



(e.g., 



A ek e et al.ll2010f ). but they are usually not seen in the spectra of T Tauri stars 



(IGeers et all 120061 ) . Therefore, for simplicity, we neglect their possible contribution to the 
opacity and heating of the disk. 

In Fig. [71 the distribution of the integrated UV field strength over the disk is shown 
i n logarithmi c scale. Contours are labeled in units of the interstellar UV radiation field 
(lDraindll978l ). In the evolved models the disk becomes more transparent to the incident UV 
radiation. While in model A5 (Fig. [7J middle panel) one can clearly see the dark zone at the 
disk midplane (up to Z/R = 0.15), this region is significantly smaller in models GS and A4. 
On the other hand, in evolved models, the midplane can be divided into two parts. The inner 
part at R < 100 AU remains as opaque to UV radiation as in the standard model, while 
the outer midplane can hardly be called "dark" , because of the smaller opacities in models 
GS and A4. In model A4, the UV strength in this region is ~ 10 5 times higher than in 
model A5. Thus, the outer disk in model A4 consists of only two layers: the atmosphere and 
the moderately UV-illuminated midplane (Fig. [TJ left panel). In the GS model the outer disk 
is not that transparent, even though the inner dark region of the midplane is still narrower 
than in model A5. The UV strength in the outer disk region is higher than in model A5 
by a factor of 10 — 10 3 , and is only a small fraction of the UV strength in the intermediate 
layer. Therefore, the three-layered disk chemical structure is preserved in the GS model but 
the midplane is only "dim" rather than "dark" (Fig. right panel). 

The UV spectrum at various heights above the midplane for a disk radius of R = 100 AU 
is shown in Fig. [8] . For illustrative purposes, we also schematically include the Lya feature. 
As expected, the dust absorption is stronger at shorter wavelengths. It is important to note 
that Lya feature "disappears" in the densest part of the disk, due to the strong absorption 
and the absence of any efficient scattering mechanisms for UV photons in our model. Hence, 
we expect only minor differences in column densities of molecules in the current version of 
our model with and without Lya, because the major contribution to column densities comes 
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from the den sest parts of the d isk. In future works, we plan to treat the scattering of Lya 
photons (see lBergin et all 120031 ). 

Another factor tha t can affect the UV field in the upper disk atmosphere is the gas- 
phase opacity. Recently, iBethell fc Berginl ( 120091 ) noted that water and OH molecules can be 
a source of significant absorption in the UV band. However, to take this effect into account 
properly, one would need to model chemical and physical disk structure self-consistently. 



4.3. Chemical structure of the disk 

We now want to discuss which effects - grain growth or sedimentation - are dominating 
changes in column densities of the various gas phase species. Note, that in all three basic 
models detailed cross-sections are used, to calculate photoreaction rates whenever possible 
(see section l3TTj) . 



4-3.1. Inner disk: models GS and A5 

First, we compare models GS and A5. The overall trend for these two models is that 
the column densities in the inner disk are greater for most species (like CO, CO2, NH3) in 
the GS model than in the A5 model, sometimes by orders of magnitude. Of all the gas- 
phase species, only 10% have greater column denisities in the A5 model at R = 10 AU (like 
S + , HNO, CP). This result may seem counterintuitive, as the disk which is transparent to 
UV radiation should represent a more hostile environment for molecules. However, at the 
same time the growth of dust particles decreases their total surface area, thus lessening the 
effectiveness of mantle formation and surface chemical processing. Also, the more intense UV 
field in the upper disk causes faster photodesorption, shifting the balance between sticking 
and desorption toward higher gas-phase abundances. 

The most dramatic example of this difference is H2CS. The column densities of this 
molecule at 10 AU in models GS and A5 differ by almost 6 orders of magnitude. How- 
ever, there are also less exotic examples represented by basic carbon compounds, also at 
10 AU, where the effects of dust evolution are stronger. Both CO (see Fig. ^ and C0 2 
are concentrated primarily in the warm molecular layer. The peak gas-phase abundances 
of CO are nearly the same in models GS and A5, as they are essentially equal to the total 
carbon abundance. However, the location of the peak of the relative abundance of CO in 
the GS model is shifted downward to the denser disk region causing the CO column density 
to increase. The peak of the CO2 abundance is also located deeper in the GS model, and 
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the maximum abundance of this molecule is two orders of magnitude higher than in the A5 
model. A similar behavior is characteristic to the abundances and column densities of N 2 H + , 
NH 3 , and H 2 (see Fig. [TD1 fourth row). While in model A5 the column densities of these 
species decrease towards the inner disk, in the models with grain growth (A4 and GS) they 
are nearly constant. Again, this is a manifestation of the faster grain evolution in the disk 
interior with respect to the outer region (R > 50 — 100 AU). 

In the inner disk midplane, where temperatures are low enough to allow icy mantle 
formation, all carbon atoms are locked in surface CO2. Midplane surface and gas-phase 
abundances of CO and CO2 do not differ in models GS and A5. At the same time, in model 
A5 surface hydrogenation products of these species are significantly more abundant. For 
example, s-H 2 CO and s-HCOOH (the prefix "s-" denotes surface species) are enhanced by 
one to two orders of magnitude in the model with "standard" dust properties. 

Another noticeable difference between models GS and A5 is the layer of enhanced de- 
pletion located right beneath the molecular layer, which is clearly observed in the A5 model, 
and is absent (or at least less pronounced) in model GS. The origin of this depletion layer is 
related to hydrogen and helium ionization. This effect is most directly seen for CO. In the 
depletion layer, this molecule is rapidly destroyed in reactions with He + . Below this layer, 
helium is not ionized as X-rays do not penetrate there and this channel is not important. 
Above it, the reaction of CO with ionized helium cannot compete with rich ion-molecular 
chemistry. In the GS model, molecular freezing is less effective because of decreased avail- 
able dust surface, and photodesorption is enhanced; thus formation of this depletion layer is 
suppressed. It is also important to note that cosmic ray (CR) and X-ray ionization is not 
affected by dust sedimentation in our model. 

Abundant gas-phase water in the GS model consumes hydrogen ions effectively (pro- 
duced by CR ionization and in reactions of H 2 with He + ), and their abundance drops signif- 
icantly at intermediate heights (about 1.5 AU at R = 10 AU). Correspondingly, abundances 
grow for other species that would otherwise have been destroyed by H + . Among these 
species are some long carbon chains, like C6H 2 , and cyanopolyynes, like HCN, HC3N and 
HC5N (Fig. (TTJ first and second rows). Conditions for the gas-phase synthesis of these com- 
pounds are only met in the narrow disk layer, and so the presence or absence of H + in this 
layer determines whether or not molecules with long carbon chains can be abundant in the 
inner disk. This is why their column densities are 4-5 orders of magnitude higher in the GS 
model, where they bind a significant fraction of C atoms, than in the A5 model. 

Thus, the typical net effect of the dust evolution is to increase the gas-phase abundances 
(despite enhanced photodissociation) of most molecules, and to decrease abundances of their 
surface counterparts. Of course, given the intricate structure of chemical networks, the 
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response is not always straightforward, and there are species for which both gas-phase and 
surface abundances are increased in the GS model (like H 2 S), and also species for which 
both gas-phase and surface abundances are decreased in the GS model (like HNO). Also, 
dust evolution suppresses surface hydrogenation and, in particular, formation of the simplest 
organic molecules in icy mantles. 

4-3.2. Outer disk: models GS and A\ 

Yet another comparison, that is useful, is to compare the results of models GS and 
A4. An increased average grain size (as well as an increased upper limit of the grain size 
distribution) can be thought of as a simplified way of modeling grain growth. The main 
assumption in this representation is that dust evolves uniformly over the entire disk. As we 
have seen, this is not the case, so the largest deviations should be in those disk regions that 
are less affected by grain growth and sedimentation, (in regions further away from the star). 
Indeed, at a radius of 10 AU, the vertical abundance profiles in the GS and A4 models for 
most species are nearly identical, with most differences observed at R > 100 AU. 

To illustrate features specific to model A4, we turn to the molecular disk content at 
550 AU. The biggest differences in column densities are found for surface species, which is 
expected, as the grain surface area is much smaller in this model than in the other two models. 
Some complex species, like s-C 2 and s-C 3 0, are underabundant in terms of column densities 
by more than 5 orders of magnitude in model A4. The column densities of other carbon- 
bearing surface species, like S-CO2, formaldehyde, methanol, are 2-4 orders of magnitude 
smaller than in other models. On the other hand, the column densities of ices, which are not 
affected that much by surface chemistry (s-CO, s-H 2 0, s-NH 3 ), differ only slightly in models 
GS and A4. This is largely because these molecules simply bind nearly all CNO atoms in 
both of these models. 

For the gas-phase species, the largest difference in column densities is observed for 
N 2 0. This molecule is overabundant (in terms of column densities) in model A4 relative 
to model GS by four orders of magnitude. Similar overabundance is also seen for other 
N-bearing species as well, like N0 2 , OCN, and ammonia. While the network of chemical 
pathways connecting all these species is quite complicated, the underlying reason for their 
higher abundance seems to be simple. Surface transformations, which are less effective at 
this radius in model A4 than in model GS, lock somewhat fewer nitrogen atoms in surface 
ammonia making them available for various gas-phase species. It must be kept in mind that 
the column density of S-NH3 is so high that even minor relative differences in iV(s-NH3) 
between models A4 and GS significantly affects column densities of other, less abundant 
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N-bearing species. Similar trends are observed in other families as well. In general, we can 
say that dust acts as an irreversible (or nearly irreversible) sink for atoms, locking them in 
surface species with high desorption energy. As grains get bigger, their total surface area 
diminishes, and more atoms are available for richer gas-phase chemistry. 

Column densities alone are not the unique indicator of a specific chemical structure. 
Even if the column densities of some molecules are almost the same in models GS and 
A4, their vertical distributions can be different in terms of both width and location of the 
molecular layer. In general, as we move from model A5 to GS to A4, at large radii, the 
molecular layer gets thinner, and is located deeper in the disk. 



4-3.3. Comparison to other works 



The chemistry in a disk with evolved dust was studied by lAikawa fc Nomural (120061 ) . In 
their investigation, Hjj" was indicated as one of a few species sensitive to the grain growth 
in terms of a colu mn density. In our study, model A4 is most suitable for the comparison 



with the model of lAikawa fc Nomural (120061 ). and comparison of models A5 and A4 indeed 
confirms the statement made in their paper. However, Hjj~ behaves differently in our two 
evolved models (GS and A4; see Fig. [101 sixth row). In model GS the column density of 
Hjj~ does not significantly differ from that in model A5. As the formation of Hjj~ is driven by 
cosmic ray ionization of H 2 , with the same rate in all three models, the difference is related 
to the destruction of this ion. It is consumed mainly in dissociative recombination reactions 
with electrons and negatively charged grains as well as in ion-molecular reactions with CO, 
H2O, CN, and other radicals. Therefore, the distribution of in the disk is controlled by 
fractional ionization and by the width and location of the molecular layer. In models GS 
and A5, peaks of molecular abundances are located at different heights, but still above the 
midplane. Because of this, the vertical profiles of the Hjj" abundance in these two models are 
similar, except for a shift relative to the midplane. Consequently, the column density of H^" 
is not sensitive to grain evolution as described by model GS. 

On the other hand, in model A4 the molecular layer extends down to the midplane, 
and, thus, protonation reactions involving H^" are active in a broader disk zone. The main 
channels of Hjj~ destruction in the outer disk are its reactions with CO and water, which 
are abundant even in the cold gas in model A4. This emphasizes the need of careful grain 
growth modeling in chemical studies, as an adopted approach to the grain evolution may 
affect conclusions even about some primary species. 



The increase in many molecular column densities seen in model GS seems to be in 
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contradiction with the conclusion of lAikawa &: Nomura! (120061 ) that grain growth does not 
affect most column densities. However, they do not consider dust settling and also model 
dust growth by increasing the upper size limit, thus, keeping the significant amount of small 
grains. This means that the upper disk in their model is, probably, more opaque than our 
disk in model GS. Thus, the increase in column densities is only observed when dust growth 
and settling are considered simultaneously. 



Our basic results compare favorably with those obtained by lJonkheid et al.l ( 12004 120071 ). 
However, the two models are also quite different. In particular, freeze-out and surface pro- 
cesses are not taken into account by Jonkheid et al., and also assumptions on the central 
star are very different from our work. 

Nevertheless, the CO depletion layer is seen in both the T Tauri disk model and es- 
pecially in the Herbig Ae/Be disk model, when dust is assumed to be well-mixed with the 
gas. As in our model, this layer disappears when dust settling is included in the model- 
ing by Jonkheid et al. In both models, the peaks in the vertical distributions of molecular 
abundances occur closer to the midplane in the models with dust evolution. 

In terms of column densities, however, the two models differ. In our model, as grains 
become larger, the mole cular column densities become larger. The opposite trend is ob- 
served in the model by iJonkheid et al.l (120071 ) (cf. their Fig. 8). But this difference is 
understandable: the disk in their model is warmer, because the central star is hotter and 
more luminous. As the disk becomes more transparent to UV radiation, molecular column 
densities drop. In our model, chemical interactions between gas and dust are very important, 
and the difference between column densities in different dust models is not controlled solely 
by photodissociation. 



Another study to be compared with our results is those of iPascucci et al.l ( 120091 ). These 
authors found that there is a significant underabundance of HCN relative to C2H2 in disks 
around cool stars and substellar objects. They attribute this underabundance to the differ- 
ences in UV radiation field. The photodissociation of molecular nitrogen is less efficient in 
the vicinity of cool stars, so there are not enough free nitrogen atoms to be incorporated into 
HCN molecules, while N atoms are abundant around hotter stars with noticeable UV excess. 
The UV spectrum of the central object is the same both in GS and A5 models, however 
one of the main results of grain evolution is that the integrated UV field in the upper disk 
is stronger in model GS than in model A5, and we may expect a similar trend in the HCN 
abundance. 



And we do see this trend in the outer disk (Fig. [T0l fifths row and Fig. [TH first row). 
While in model A5 at R > 100 AU column densities of HCN and C2H2 are almost equal, in 
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model GS (with stronger UV field) HCN is 15 times more abundant than acetylene. Analysis 
of chemical reactions shows that enhanced HCN is indeed caused by greater abundance of 
N atoms. The trend is not preserved at smaller radii, though. At R = 10 AU in model 
A5 the column density of HCN exceeds that of acetylene by a factor of 42 (compared to 
near equality at greater distances). This is because in model A5, with abundant dust in the 
upper disk, surface HCN synthesis with subsequent desorption becomes an important factor 
increasing gas-phase abundance of this molecule. 

One should keep in mind that both under- and overabundance of HCN relative to C2H2 
in model GS is observed on top of the overall trend for column densities of most molecules 
to be greater in this model. 

4-3.4- Fractional ionization in the disk 

Dust evolution is reflected not only in the abundances of some molecules, but also 
shows up in a more integral way, affecting the fractional ionization x e of the disk (see Fig. 
[TOj bottom row). Since in our model one grain can adsorb not more than one electron, change 
of dust properties impacts disk ionization mainly through change of absorption of the UV 
field which is important for gas-phase chemistry. While the ionization degree exceeds 10~ 10 
everywhere in the disk, grain number density is not more than 10 -12 in respect to hydrogen 
nuclei. Therefore, ionization degree is controlled by cation chemistry while the contribution 
of electron sticking to dust grains is less than 1%. In the disk atmosphere differences are 
small as the entire ion content is provided by ionized carbon and ionized hydrogen (in the 
uppermost X-ray ionized layer). Deeper in the disk, different properties of dust cause x e in 
models A5 and GS to differ by almost two orders of magnitude. This difference is primarily 
caused by different dominant ions at different heights. Let's consider a slice at R = 100 AU 
as an example. When absorption gets high enough so that C + is no more a dominant ion 
(z < 40 AU), in model GS most electrons are provided by metals. Sodium, magnesium, 
and iron ions stay dominant down to z ~ 15 AU. Below this level, where most dust is 
concentrated in model GS, metals are depleted, and most abundant ions are H + and H^" (in 
almost equal shares). 

In model A5 dust is equally abundant everywhere in the disk, and metals are depleted 
at all heights below 40 AU. Correspondingly, almost all ions at these heights are provided 
by H + and Hj. These species are heavily involved in ion- molecular chemistry, so their abun- 
dances are very sensitive to the overall molecular content. In particular, ionized hydrogen 
is effectively consumed in reactions with ammonia, water, hydroxyl etc. These species are 
abundant in the warm molecular layer (z ~ 30 AU) and depleted right below it (z ~ 20 AU). 



-24- 



Correspondingly, in the warm molecular layer Hg is the most abundant ion, while at lower 
heights the balance shifts toward H + . The equilibrium fractional ionization (for a single ion) 
is inversely proportional to the square root of the recombination coefficient. This coefficient 
is significantly lower for H + than for ITj", so the fractional ionization is higher at z ~ 20 AU, 
where H + dominates. 



Discussion 



5.1. Disk chemistry as a tracer of grain evolution 



Can molecular abundances in a protoplanetary disk be tracers of grain evolution? Since 
the distance to the Taurus-Auriga complex, where mos t known disks around low-mass pre- 
main sequence stars are located, is ~ 140 pc (lEliaslll978l ) and the expected angular resolution 
of ALMA is up to 0".01, the spatial resolution of disk observations should be ~ 1-2 AU. 
This is sufficient to study the spatial distribution of molecular abundances in the outer disk 
(10-600 AU). So it is worthwhile to ask if there are species, sensitive to the grain growth at 
a level that would make them detectable tracers of this process. An immediate product of 
observations are molecular lines. Therefore, we need to predict not abundances or column 
densities, but the shapes and intensities of m olecular lines, as well as the variations caused 
by grain evolution (e.g., ISemenov et al.ll2008l ). The physical model of the disk employed in 



this study is quite simple, however, so it does not make sense to perform realistic radiative 
transfer modeling. In this situation, we have chosen to consider only given species as a 
possible grain evolution tracer if its column density differs in models A5 and GS by an order 
of magnitude or more in at least one of representative radial regions around 10 AU, 100 AU 
and 550 AU. An exception is made for CO, which has a very high column density. In the 
inner disk and at R = 100 AU from the star CO column density becomes higher by a factor 
of 5-7 in model GS in comparison to model A5. 

Major species, both sensitive and insensitive to grain evolution in our model, are grouped 
in Tables [2] and [3] respectively. Additionally, Table [2] shows the observed column densities 



for several species, taken from ISemenov et a 



(120041 ) (except for HCO+ an d N^H+ which 



are ta ken from Table 3 in iDutrey et al.l (120071 ) and C2H which is taken from iHenning et al. 



(120 101 )). It is clear that the column densities of species which have already been observed 
in disks so far, do not seem to be strongly influenced by the changing dust properties at 
large distances from the star, i.e., in the domain of single-dish observations. Closer to the 
star, however, the differences are more significant. Specifically, at R=10 AU, the column 
density of water in model GS exceeds that in model A5 by a factor of 560. A word of caution 
is to be said about water abundance. While column densities for all other molecules have 



-25 - 



reasonable values in our models (as co mpared to availabl e observations), water seems to be 
an exception. Recent H erschel results (|Bergin et al.ll2010f) show no observable water lines in 
the DM Tau spectrum. iBergin et al.l (120101 ) estimated the upper limit for the water column 



density in the outer disk to be about 3 ■ 10 13 cm" 2 , which is much lower than the water 
column density in our models. They suggested that low water abundance may serve as an 
indication of grain growth and sedimentation. However, as we said previously, our model do 
not predict significant grain evolution in the outer disk, so the DM Tau disk may represent 
a later stage of the grain growth than our model. 

Of course, not only abundances (i.e., column densities) themselves, but also their ratios, 
can be an indicator of the grain growth. However, as in the example of HCN and C2H2, 
considered above, and motivated by observations by lPascucci et al.l ( 120091 ). the influence of 
the grain evolution on the abundance ratio can be at least twofold, and calls for a more 
detailed consideration. 

Another pair of observed s pecies that may pot entially be useful as a tracer of the grain 
evolution is HCO + and N 2 H + ( iDutrey et al.l 120071 ). While the av erage column densit ies of 
these species in model GS agree quite well with observations of IDutrey et al.l ( 120071 ). the 
ratios of their column densities in models GS and A4 are quite different: at R = 550 AU the 
ratio A^(HCO + )/A^(N 2 H + ) is 150 in model GS, and 1500 in model A4. While the main focus 
of our work was on comparison between models GS and A5, not between models GS and 
A4, the difference mentioned here shows that this ratio is worth checking in a more detailed 
model. 

At the disk periphery, where dust evolution is least prominent, the only species that is 
sensitive to dust growth and sedimentation on the level of an order of magnitude is methanol. 
However, the methanol column density is quite modest (2-10 11 cm -2 ) and even smaller closer 
to the star. 

In general, we may conclude that the chemical signatures of early grain evolution are 
hardly observable with single-dish instruments, although they are a promising target for 
highly sensitive interferometers such as ALMA. 

The composition of the icy mantles of dust grains in the d isk also may provide infor- 
mation about the grain evolution process. iTerada et al.l (120071 ) showed that it is possible 
to assess the mantle composition through observations of scattered light in edge-on disks 
around T Tau stars. They found deep water-ice absorption lines corresponding to an optical 
depth of r=l in the directions of edge-on disks HK Tau B and HV Tau C. In order to check 
whether ices in the observed region are sensitive to the grain evolution process, we calculated 
the mantle composition in models A5 and GS near the optical disk surface r = 1 (see Table 
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H]). Due to the different dust properties in models A5 and GS, the location of this surface 
is different in the two models. It is closer to the midplane in model GS, where the disk is 
more transparent to UV radiation than in model A5, where the dust is pristine and more 
opaque. The vertical temperature gradient of the dust is the most important factor in de- 
termining the differences in the ice composition between the two models. Since our model of 
the vertical structure of the disk is not self-consistent, the temperature at the r = 1 surface 
is determined only approximately. Therefore, differences in the chemical composition of ices 
in models A5 and GS should only be considered as qualitative in nature. 

In the inner disk (R = 10 AU), the temperature at its optical surface is 50 K in model GS, 
and 80 K in model A5. At this temperature range, the ice composition is totally dominated 
by water (~ 99%). Ammonia is the second most abundant ice compound (~ 1%), and 
other species are present in ice only in trace amounts. At a disk radius of 100 AU, the ice 
composition at r = 1 is significantly different in models A5 and GS and relatively complex. 
While water is still the dominant ice constituent, it shares this role with hydrocarbons C n H m , 
which are produced efficiently by both surface chemistry and gas-phase chemistry in this disk 
region. Also, in the GS model, several other species, such as C0 2 , H 2 CO, HCN etc. stay 
in ice in noticeable amounts. Further away from the star, at a radius of 550 AU, the ice 
composition in models A5 and GS is almost the same again because the optical surface at this 
distance in both models corresponds to very similar temperatures of approximately 30 K. 

Such differences in grain mantle composition at the optical surface of the disk may serve 
as an additional tracer of the grain evolution process. However, the results described above 
should be confirmed with detailed, self-consistent models of the vertical structure of the disk. 

5.2. Model limitations and future work 

We have already mentioned that column densities are only a partial indicator of grain 
growth. Line emission modeling is needed to assess observable tracers of dust evolution 
with more confidence. But, there are some issues that need to be addressed prior to this 
modeling. As we mentioned in the introduction, the simplified setup of our model allows 
us to isolate the chemical outcomes of grain evolution, although the process of grain growth 
should probably have some more general influence on the disk structure. The response of 
a disk to dust evolution has already been addressed a number of times but with only an 
oversimplified approach to grain growth and/or sedimentation. Our modeling shows that 
dust evolution in a disk is highly non-uniform, both in the vertical and radial direction, 
which should be reflected in disk physical properties. 
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In future disk models, one would need to at least calculate the processes of grain growth 
and sedimentation consistently with the description of the disk structure, as density and 
temperature distributions in the disk are controlled by opacities. Dust absorption is definitely 
affected by growth and sedimentation, but both these processes, in turn, depend on the disk 
structure. In this study, changes in the surface- averaged grain size (and, hence, in the 
total grain surface area) are only minor, and even in model GS, where dust is evolved, 
the disk in relatively opaque to the incoming radiation. This means that extra heating 
due to decreased opacity would mostly affect the upper disk, above the region where most 
molecules are concentrated. As in this region molecular abundances are mostly controlled 
by photoprocesses, most column densities probably will not be strongly affected by the 
increased gas scale height or gas temperature in the upper disk. The exception is represented 
by abundant components that bind almost all CNO atoms in the molecular layer (CO, N 2 , 
O). As we noted, the increase of their column densities in the inner disk in models GS and 
A4 relative to A5 is partially caused by the downward shift of the molecular layer to the 
denser disk region. If we would take the disk vertical expansion due to increased heating 
into account, the peak of relative abundance for these components would probably reside in 
the region of smaller density, thus, suppressing the increase in column densities. This calls 
into question the ability of CO to be a tracer of dust evolution within the fram ework of our 



mode l. However, dust growth may affect not CO itself, but its isotopic ratios ( jVisser et al. 



200% . 



Another issue, which should be addressed properly in the future, is the possible in- 
terrelation between chemical and physical structure of the disk with evolving dust. This 
interrelati on has already been discussed from the point of view of thermal balance in the 



disk (e.g.. |Jonkheid et al.ll2007l ). However, dust evolution may also lead to significant changes 



in the fractional ionization, which is an important dynamical par ameter, especially if tu rbu 



lence in disks is excited due to the magnetorotational instability ( lllgner fc Nelson! 120081 ) . All 
of this calls for more self-consistency not only between disk structure and grain evolution, 
but also between these two factors and chemistry. 

The next step is to include time dependence in the model. Two of the three ingredients 
of our model (disk structure and grain ensemble properties) explicitly imply that the disk 
is in a steady state. The third ingredient (chemistry) is formally time-dependent, but pro- 
duces steady-state abundances for most species. This is the reasonable outcome, given the 
equilibrium nature of the other two ingredients. Our results represent only a snapshot of the 
disk at an early (or even the earliest) stage of grain evolution. However, at some later stage, 
the grain growth timescale may become short enough to be comparable with the timescale 
of various chemical transformations. With the current model we cannot say what impact 
on chemistry possible further grain growth, beyond what the Birnstiel model can currently 
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calculate, would have, even though some hints are provided by model A4. In particular, 
we may expect higher molecular densities at the disk periphery, with the molecular layer 
extending down to the midplane. This confirms the earlier suggestion that observations of 
gas-phase CO and water molecules in a cold disk region can be explained by an increase 
in grain size. Further improvements in the dust growth model are also needed. This field 
now represents an area of active research, both theoretical and experimental, and there are 
certa in indications tha t grain evolution is not limited to simply sticking and fragmenting 



ve-g- 



Zsom et al.ll2010h . 



Another issue that may become more important with the grain ensemble containing 
large grains is the treatment of grain surface chemistry. Currently, in chemical modeling we 
consider grains of equal "effective mean size" instead of the real grain size distribution. We 
believe that this simplification does not dramatically affect results of this study. The average 
grain size defined according to (Q is close to the size of grains that dominate the total grain 
surface area available for surface chemistry. In other words, most of the surface chemistry 
occurs on grains of similar sizes. Also, Acharyya & Herbst (in preparation) constructed 
a chemical model with grains of several grain sizes in order to reproduce the grain size 
distribution. They found its results to be quite similar to those obtained with a single grain 
size. However, the problem may be more complicated as the temperature of an individual 
grain depends on its size. For smallest grains, stochastic heating becomes important, which 
makes it impossible to assign a single temperature to them. These questions are not trivial 
and worth to be considered in a separate study. In this work, we do not have grains of 
extremely small radii (R < 10~ 6 cm). So, the considerations above should be of only a 
limited importance here. 

Finally, transport processes are an important evolutionary factor in an accretion disk. 
At least, some of the differences between models GS and A5, like the enhanced depletion 
layer in model A5, are quite localized. One does need to account for possible gas mixing to 
check whether these differences can survive in a more dynamic medium. The radial motion 
of grains has to be included in a realistic treatment of the disk evolution. 



6. Conclusions 



We investigated the chemical evolution of the protoplanetary disk around a T Tauri 
star. The thermal and density structures of the disk are calculated with a 1+1D a-viscosity 
model. The pro cesses of grain coagul ation, fragmentation, and se dimentation were modeled 
as described in iBrauer et al.l (l2008bl ) and iBirnstiel et al.l (120101) . T he initial dust ensem- 



ble is the standard MRN grain size distribution (IMathis et al. 



19771 ). Using the final grain 
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size distribution (at t — 2 Myrs) and the 1+1D disk structure, we calculated the verti- 
cal distribution of the dust in the disk, assuming stirring-settling equilibrium according to 



Dullemond &: Dominikl ( 120041 ) . For this distribution, we calculated the approximate UV field 
in the disk, and simulated the chemical evolution of the disk over 2 Myr. Results for three 
dust models were compared: (1) a classical dust model, with uniform spherical particles of 
0.1//m and a gas-to-dust mass ratio of 100 (model A5); (2) a simplified grain growth model, 
with larger grains of 1/xm (model A4); (3) a detailed grain growth model (model GS). The 
main results of the work can be summarized as follows: 

• The fragmentation of grains due to collisions keeps a population of small dust particles 
in the disk, which dominate the total grain surface area crucial for chemical evolution 
of the disk. This effect reduces the impact of grain growth on the chemical structure 
of the disk. Grain growth proceeds fast in low-mass T Tauri disks at distances < 
20 — 50 AU where the steady-state grain size distribution is reached in a few thousand 
years. Outside this radius, up to 10 6 years is needed to reach steady state at R>100 AU. 
Furthermore, grain coagulation is limited there, and increases the average grain size 
by no more than an order of magnitude in the midplane. After 2 Myrs of evolution, 
the average grain size in the disk midplane varies between 10~ 4 cm at R = 10 AU, and 
2 • 10~ 5 cm at R = 550 AU. The dust-to-gas mass ratios in the same distance range 
are confined between ~ 10~ 6 and 0.02. 

• The net effect of the dust settling and growth is a reduction of the total grain surface, 
and higher UV irradiation rates in the upper disk. The chemical structure of an 
evolved disk still has three layers, but the intermediate molecular layer gets wider, and 
shifts closer to the midplane. Therefore, the abundances and column densities of many 
species are enhanced by a factor of 3 — 100, even for such a moderate grain growth. 

• A simplified model of grain growth, in which the dust-to-gas mass ratio is kept constant 
and grain size is simply increased by one order of magnitude, is not sufficient to re- 
produce the chemical evolution of a disk with evolving dust. For e xample, the column 
densi ty of H^, which is proposed to be sensitive to grain growth (lAikawa &: Nomura 



20061 ) . exhibits a high sensitivity to grain growth in our simple model, but not in the 



model with a detailed treatment of grain growth and sedimentation. 

• Grain evolution suppresses the formation of organic molecules in icy mantles, but also 
enhances gas-phase molecular abundances. In particular, water and CO abundances 
are enchanced in the cold midplane. 



• We propose a few observational tracers of grain evolution process. These are the column 
densities of gas-phase molecules, such as C 2 H, HC 2n +iN (n — — 3), H 2 for the inner 
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disk. The abundance ratio of C 2 H 2 /HCN may serve as a grain evolution tracer in 
outer disk, as well as the abundance of CH 3 OH. The composition of the ice mantles of 
interstellar grains seen in scattered light in the intermediate disk layer (around r=l) 
may also serve as a tracer of grain evolution. 
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Table 1: Main disk model parameters 



Inner radius 


0.03 AU 


Outer radius 


550 AU 


a 


0.01 


M, M year 1 


2 • 10~ 9 


X-ray luminosity 


10 30 erg s- 1 


Cosmic ray ionization rate 


1.3 • 10- 17 s- 1 


Disk mass 


0.055 M 


Initial average grain size 


10~ 5 cm 



Table 2 : Species sens i tive t o grain evolutio n. Ob served column dens ities are compiled from iDutrev et al.l ( 119971 ); |Q . 

J2000h : lAikawa et all J2OO2I ): butrev et all J2007h . bergin et all fcoioh and lHenning et all J201ph . 



Species 






Column densities, cm 


-2 






Peak abundance 


Observed 




10 AU 




100 AU 


550 AU 




n(X)/n(H) 


column densities, cm -2 




A5 


GS 


A5 


GS 


A5 


GS 


A5 


GS 


DM Tau 


LkCal5 


CO 


2.0(17) 


1.1 


18) 


2.0(17) 


9.4(17) 


1.7(17) 


2.9 


;i7) 


7.3(-05) 


7.3(-05) 


5.7(16) 


9.0(17) 


co 2 


6.8(10) 


9.5 


13) 


2.4(12) 


8.5(14) 


8.2(13) 


9.8 


:i4) 


7.2(-08) 


4.6(-07) 


- 


- 


CH30H 


1.1(03) 


3.5 


;os) 


1.9(08) 


1.0(10) 


4.6(09) 


1.4 


10) 


1.2(-11) 


1.2(-11) 


- 


- 


H 2 


2.5(13) 


1.1 


16) 


3.6(14) 


7.0(15) 


1.4(14) 


2.0 


;i5) 


1.3(-06) 


2.6(-06) 


<3.0(13) 


- 


H 2 S 


1.7(05) 


1.3 


10) 


2.0(09) 


1.1(11) 


2.9(10) 


3.6 


10) 


2.0(-ll) 


1.9(-11) 


- 


- 


C 2 H 


6.3(10) 


2.0 


12) 


2.6(12) 


2.2(12) 


6.7(12) 


5.0 


12) 


1.4(-09) 


1.5(-09) 


2.8(13) 


2.9(13) 


C 2 H 2 


6.2(10) 


4.6 


12) 


4.5(12) 


3.0(12) 


4.4(12) 


1.6 


12) 


4.0(-09) 


1.3(-09) 


- 


- 


CH3CH3 


3.0(09) 


1.5 


12) 


1.8(10) 


3.3(08) 


1.3(07) 


2.5 


;o6) 


3.6(-ll) 


2.3(-10) 


- 


- 


H 2 CS 


8.0(05) 


3.5 


11) 


2.7(11) 


8.3(11) 


2.7(10) 


1.8 


;n) 


2.2(-10) 


5.6(-10) 


- 


- 


HON 


2.6(12) 


5.0 


13) 


6.9(12) 


4.2(13) 


9.0(12) 


2.1 


13) 


1.8(-08) 


9.1(-09) 


2.1(12) 


7.8(13) 


HC 3 N 


6.8(08) 


1.7 


12) 


4.1(11) 


2.8(11) 


1.1(12) 


4.3 


10) 


9.1(-10) 


l.l(-09) 






HC 5 N 


4.9(08) 


1.5 


12) 


2.1(11) 


1.5(11) 


1.8(11) 


8.8 


10) 


5.9(-10) 


6.5(-10) 






HC 7 N 


4.7(06) 


1.2 


12) 


6.1(10) 


5.3(10) 


9.2(10) 


2.8 


10) 


3.9(-10) 


6.0(-01) 






HCNH+ 


2.0(10) 


2.7 


11) 


5.7(10) 


1.3(11) 


6.6(10) 


6.8 


10) 


1.1(-10) 


3.6(-ll) 






HCOOH 


1.4(11) 


6.3 


13) 


1.1(12) 


1.3(13) 


2.6(11) 


2.1 


12) 


4.0(-09) 


1.2(-08) 






OCN 


3.2(07) 


4.3 


;o9) 


2.3(10) 


1.5(12) 


1.3(11) 


4.2 


12) 


2.8(-10) 


4.3(-09) 






OCS 


1.6(07) 


9.3 


10) 


1.6(10) 


1.5(10) 


2.7(08) 


4.5 


;o8) 


1.9(-11) 


1.9(-11) 




<2.9(13) 


NH 3 


1.7(11) 


3.2 


13) 


1.1(13) 


3.5(13) 


8.5(12) 


1.2 


13) 


1.3(-08) 


9.4(-09) 






HCO+ 


1.1(11) 


5.2 


12) 


2.6(12) 


2.2(12) 


2.2(12) 


1.9 


12) 


2.0(-09) 


7.0(-10) 


6.5(12) 


8.0(12) 


OH 


6.1(13) 


1.9 


14) 


2.7(13) 


1.2(14) 


1.9(13) 


1.5 


;i4) 


3.4(-08) 


7.0(-08) 







(2000 


) and 


Dutrev et al. 


(2007) 



Species 






Column densities 






Peak abundance 


Observed 




10 AU 


100 AU 


550 AU 






column densities, cm -2 




A5 


GS 


A5 GS 


A5 


GS 


A5 


GS 


DM Tau 


LkCal5 


CN 


6.8(12) 


2.9(13) 


5.6(12) 9.6(12) 


1.2(13) 


1.7(13) 


3.6(-09) 


2.7(-09) 


9.5-12(12) 


6.3(14) 


C 


5.1(16) 


5.9(15) 


6.3(16) 7.4(16) 


1.6(17) 


2.1(17) 


5.6(-05) 


5.4(-05) 






H 3 + 


2.3(12) 


3.7(12) 


8.5(12) 1.2(13) 


2.0(13) 


2.5(13) 


2.4(-09) 


1.0(-09) 






N 2 H+ 


3.0(09) 


2.4(10) 


1.7(10) 1.7(10) 


1.9(10) 


2.2(10) 


6.1(-11) 


2.4(-ll) 


1.1(11) 


2.3(11) 


HNC 


1.9(12) 


4.3(12) 


2.5(12) 1.5(13) 


5.4(12) 


1.3(13) 


6.3(-09) 


4.7(-09) 


9.1(11) 


<5.4(12) 


C+ 


2.5(16) 


6.2(16) 


3.3(16) 7.8(16) 


5.0(16) 


8.2(16) 


7.3(-05) 


7.3(-05) 






CH 3 


2.8(11) 


1.7(12) 


4.7(12) 3.0(12) 


2.4(13) 


1.8(13) 


7.0(-09) 


1.7(-09) 
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Table 4: Ice composition in the disk at the surface r = 1. Only species with contribution of 
1% or more are shown 



Species 


Ice fraction, % 








10 AU 100 AU 


550 AU 




A5 GS A5 GS 


A5 


GS 


s-H 2 


99.8 97.3 57.4 91.0 


70.0 


88.9 


s-NH 3 


1.0 - 


1.0 


5.0 


s-C 2 H 2 




4.6 


1.0 


s-C 3 H 2 


20.9 3.3 


16.5 


1.5 


s-C 5 H 2 


12.0 - 


2.1 




s-C 7 H 2 


5.0 - 






s-C 9 H 2 


1.3 - 






s-C0 2 


2.0 






s-H 2 2 






1.5 


s-HCN 




3.1 


1.0 


s-HNO 


1.0 






s-H 2 C 4 


1.0 - 


1.0 
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Fig. 1. — The solid line shows the radial dependence of the gas surface density in the adopted 
disk model. The dashed line represents a power-law fit with the slope of -0.85. 
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Fig. 3. — Evolution of surface density of growing dust in model GS. Contours show surface 
density in gem -2 for dust grains of a certain size at a certain radius, so that, e.g., surface 
density of 0.1 fim grains at 1 AU is about 10~ 3 gem -2 after 10 4 years of evolution (leftmost 
panel). Dust growth is faster in the dense inner part of the disk. Steady-state grain size 
distribution is reached there between 10 4 and 10 5 years. In the outer disk, density is lower 
and steady-state distribution is reached between 10 6 and 2-10 6 years. 
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Fig. 4. — Radial distribution of the vertically-integrated average grain radius at different 
time moments. Averaging is performed in a way that preserves the total grain mass and 
surface area. 
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Fig. 5. — Dust-to-gas mass ratio (upper row) and an average grain radius (bottom row) in 
model of the disk with grain growth and sedimentation. At the left panels, 2D distributions 
over the disk are shown while at the right panels vertical distributions at three representative 
radii are plotted. Large grains that contain most of dust mass sediment to the midplane 
more efficiently than small ones. This reduces the dust-to-gas mass ratio above the midplane 
and decreases the average grain size in upper layers of the disk. 
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GS / A5 grain surface area ratio 




10 100 
R, AU 

Fig. 6. — Ratio of grain surface areas in models with- and without grain evolution. Grain 
evolution reduces the total grain surface area. This effect is more significant in the inner 
disk because grain growth is more efficient there, and also above the midplane because of 
grain sedimentation to the midplane. 
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10 100 10 100 10 100 

R, AU 

Fig. 7. — Logarithm of the UV field intensity in the disk for models A5, A4 and GS measured 
in units of the interstellar UV radiation field. Changes in grain properties strongly affect 
the amount of UV radiation penetrating into the disk. While A5 model has a dark midplane 
area at all radii, the model A4 disk is relatively transparent to the UV radiation beyond 
100 AU. In the model GS the extent of the dark region is smaller than in model A5 (for 
R<100 AU) and there is a "dim" rather than a "dark" midplane beyond 100 AU. 
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Fig. 8. — Spectrum of the UV radiation in the disk at different heights above the midplane 
at R = 100 AU in model GS. Closer to the midplane, the spectrum shape is similar to that 
of the ISRF because the stellar component is attenuated faster than the interstellar one. For 
comparison, spectrum of the scaled-up Draine field is also shown. 
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Fig. 10. — Vertical distributions at 10, 100, and 550 AU and column densities of CO, HCO + , 
C 2 H, H 2 0, C 2 H 2 , H3 , and electrons in models A5 (dashed line), GS (solid line), and A4 
(dotted line). 



-50 - 




Fig. 11.— Same as in Fig. EH but for HCN, HC 3 N, CH 3 OH, C0 2 , s-NH 3 , s-C 3 H 2 , s-H 2 0. 



